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ABSTRACT 

We present an analytical approach to calculate the hydrodynamics of the interaction between 
a relativistic ejecta and a surrounding medium, whose evolution serves as a model for Gamma- Ray 
Burst afterglows. We investigate the effect of the relevant model parameters on the X-ray, optical 
and radio fluxes, and the effect of a refreshed shock energy input and anisotropy in the ejecta 
on the shape of the light-curves. We compare our numerical results to observed afterglows and 
give a quantitative description of the conditions (geometry and physical parameters) in the ejecta 
that are compatible with the light-curves of the 970508 afterglow, for which a large number of 
accurate flux measurements are available. We find that the radio, optical and X-ray light-curves 
of this afterglow can be explained satisfactorily within the spherically symmetric fireball model, 
assuming a delayed energy injection, or by an axially symmetric jet surrounded by a less energetic 
outflow. 

Subject headings: gamma-rays: bursts - methods: numerical - radiation mechanisms: non-thermal 

1. Introduction 

Afterglows from Gamma-Ray Bursts (GRBs) have been observed from a number of objects at X-ray, 
optical, and in one case also at radio wavelengths. Simple analytical models are successful at explaining the 
major features of the light-curves (Meszaros & Rees 1997, Vietri 1997, Tavani 1997, Waxman 1997; Wijers, 
Rees & Meszaros 1997). The optical and X-ray light curves presented by many authors (e.g. Pedersen et 
al. 1998, Piro et al. 1998, Garcia et al. 1998, Bartolini et al. 1998) have provided evidence for occasional 
departures from the basic overall power-law decay behavior. Such departures, as well as the possibility 
of temporal power-law decays that are not exclusively determined by the spectral index, have been shown 
to follow naturally from fireball models where the radiative regime changes, the energy is not distributed 
isotropically in the ejecta (Meszaros , Rees & Wijers 1998), or where the energy input depends on the 
Lorentz factor during the brief injection episode of the central engine, leading to refreshed shocks (Rees & 
Meszaros 1998). Here we go beyond simple analytical asymptotic models, we derive and solve numerically 
the differential equations for the dynamics of the afterglow in the general case of a inhomogeneous external 
medium and refreshed shock mechanism, and calculate numerically the light-curves arising in such scenarios. 

Our previous numerical work (Panaitescu & Meszaros 1997,1998) on simulations of light-curves and 
spectra was based on a hydrodynamic code (Wen, Panaitescu & Laguna 1997) that solves the equations of 
relativistic hydrodynamics and and the shock jump conditions. The energy release mechanisms (synchrotron 



-2 - 



and inverse Compton) were treated as described in Panaitescu & Meszaros (1997). A calculation of the 
spectra and time history of an afterglow from a spherically symmetric shocked fireball is equivalent to 
computing a quadruple integral: over the lab frame time, over the structure of the shocked fluid, over 
the angle relative to the line of sight toward the fireball center (LSC) of symmetry and over the electron 
distribution. The hydrodynamic timesteps required to propagate the shell of shocked fluid over times that 
are more than 5 orders of magnitude larger than the shell crossing time and those necessary for an accurate 
calculation of the radiative losses lead to exceedingly long numerical runs, which are not best suited for an 
investigation of the effects of the large number of model parameters involved in the typical external shock 
scenario of GRBs and afterglows (Meszaros & Rees 1997). The numerical task is even more time-consuming 
in the case of anisotropic ejecta, where a new integral over the azimuthal angle is added. 

To acquire computational speed, we have developed a numerical code that calculates accurately the 
evolution of the remnant shell's flow Lorentz factor, by solving the equation that gives the evolution of the 
kinetic energy of the remnant during the ejecta-external medium interaction, with allowance for an energy 
injection in the reverse shock (Rees & Meszaros 1998), and radiative and adiabatic losses. Anisotropy of 
the ejecta or of the energy input is included at the simplest level, assuming cylindrical symmetry around an 
axis that is not necessarily the same as the LSC. Despite the assumed degree of symmetry in the ejecta, the 
resulting light-curves show a great diversity. Possible inhomogeneity of the external medium is considered in 
the form of a power-law density. To simplify the energy release treatment, we ignore here the inverse Compton 
scattering of the self-generated synchrotron photons, which is a fairly good assumption, substantiated by 
our previous results (Panaitescu & Meszaros 1998). In what follows we describe the analytic treatment of 
the remnant's dynamics and energetics, derive analytic light-curves and present our numerical results. We 
discuss the effect of model parameters on the features of the numerical light-curves, and compare them with 
the afterglow of GRB 970508. 

2. Hydrodynamics of the Remnant and Energy Release 

The most important parameter characterizing the temporal evolution of the afterglow is the bulk Lorentz 
factor r of the contact discontinuity between the ejecta and the swept up external matter. The evolution of 
L is determined by two main factors, the hydrodynamics of the shell (including the energy input, adiabatic 
losses and the deceleration caused by the external medium) and the radiative losses (synchrotron cooling). 

2.1. Adiabatic remnant 

In the absence of a delayed energy injection and of radiative losses (adiabatic remnant), at any time, 
the total energy of the fireball is constant: d[M(L — 1) + TU] = 0, where M = M + M ex is the total remnant 
mass (the sum of the initial ejecta mass Mo and the swept-up mass M ex ) and U is the co- moving frame 
internal energy of the remnant. The evolution of U is given by the adiabatic losses and by the heating of the 
external matter: dU = (dU) a d + (dU) ex . The jump conditions (Blandford & McKee 1976) at the forward 
shock imply that (dU) ex = (T — 1) dM ex c 2 , where c is the speed of light, therefore the energy conservation 
can be written as: 

(Mc 2 + U) dT + T(dU) ad + (L 2 - 1) dM e ,c 2 - , (1) 

where dM ex = 4irr 2 p ex (r) dr and (dU) a d = —(7 — 1)U d(\nV') a d- The first term in equation (|l]) is the 
lab-frame change in the kinetic energy, the second represents the adiabatic losses and the third is the total 
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lab- frame kinetic energy of the shocked external medium, including its internal energy. In the above equations 
r is the radial coordinate of the thin remnant, p ex — Pd{r /rd)~ a , Pd is the external medium density at the 
deceleration radius r^, defined as the radius where the swept up mass is a fraction T^" 1 of the initial fireball 
mass M , T being its initial Lorentz factor, a < 3 is the index of the external matter power-law density, 
7 is the adiabatic index (maintained close to 4/3 even when T < 2 by the relativistic electrons, provided 
that the fraction of electrons that are shock-accelerated is not too much below unity) and V is the remnant 
comoving volume. The differential d(ln V') a£ ; refers to that part of the comoving volume that is occupied by 
the already shocked fluid, i.e. it excludes the change in the comoving volume due to the sweeping up of the 
infinitesimal dM ex . 

We assume that the shocked external fluid stores most of the internal energy of the remnant and that 
it gives most of the afterglow radiation, as the forward shock is more relativistic than the quasi-newtonian 
reverse shock, and thus more efficient in converting kinetic energy into internal and in accelerating high 
energy electrons that are less adiabatic than the electrons accelerated by the reverse shock. Therefore, in 
this work we neglect the dynamical and radiative importance of the reverse shock, and leave its treatment for 
a future, more detailed study. We also assume that the remnant volume is practically given by the volume of 
the shocked external medium (which is correct if the injected mass is not too large compared to the mass of 
swept up external medium, because the ejecta compressed between the contact discontinuity and the reverse 
shock is denser than the shocked external medium) and that the comoving density behind the blast wave is 
equal to that set by the shock jump conditions. This implies that the comoving density is determined solely 
by p ex (r) and T: p'(r) = (7 + + 1) p ex {r). Using d(\nV') a d = — d(ln p'), the adiabatic losses can be 

written as: 

(dto— -tt-D^-^f)^. (2) 

The shocked external medium mass is given by 

dM ex = 3-^— dr, (3) 
1 r d 

which allows one to calculate the comoving volume V = M ex /p'. 



2.2. Delayed Energy Input 

It is possible that the material injected by the cataclysmic event that generates the relativistic fireball 
does not have a unique To, and that some material is ejected with lower initial bulk Lorentz factors, down to 
some limiting T m . Following Rees & Meszaros (1998), we shall consider that all the ejecta has been released 
impulsively (on a timescalc short compared to the afterglow timcscale), all at the same location, and with 
a power-law distribution of energy per unit Lorentz factor F^-: dEf cx r^Mr^, for F m < Tf < Yd, where 
Td < To is the Lorentz factor of the contact discontinuity at r = rd (for reasons given below, we start our 
simulations at r = rd). The constant of proportionality is determined by the total injected energy E in j, 
which will be one of the free parameters of the model. The fluid moving at lower T f lags behind the contact 
discontinuity and catches up with it later, as the fireball is progressively decelerated by the interaction with 
the external fluid. From the kinematics of the problem, the Lorentz factor of the ejecta that interacts with 
the reverse shock at radius r is given by 
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(r,/r) 2 -i ^. (4) 
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The mass injected in the remnant up to radius r satisfies the differential equation 



dr 

dM f =F s (T f ,T)M inj — , (5) 



where 

Fir, n = 

2 i-(r m /r d )* Vr/ 

Here Mi n j is the total mass that is eventually injected in the remnant, which corresponds to an Ei n j given 

by 

s(s - 1)- 1 [i - ( Tm/Tdy- 1 ] r m + (r m /r d y 
i - {r m /r d y 



■E inJ - = fc J M mJ C . (7) 



As the shocked shell propagates from r to r + dr, the infinitesimal injected mass dMf given by equation 
(^) collides with the shell, increasing the remnant kinetic energy by dE^j and its internal energy by (dll)f. 
These infinitesimal energies can be determined from momentum and energy conservation: 

dE kJ = r f [i-T 2 (i-/3p f )]dM f <? , (8) 
(df/) / = [rr / (i-/? / 3 / )-i]dAf /C 2 , (9) 

where [3 represents the velocity in units of c. 



2.3. Radiative Losses 



As mentioned before, we consider here that the shocked fluid cools through adiabatic expansion and 
emission of synchrotron radiation from electrons accelerated by the forward shock. There could be some 
contribution to the early afterglow light-curve from electrons accelerated by the reverse shock, but this is soon 
overcome by the forward shock emission that shifts toward lower energies, as the remnant is decelerated. 
We assume that nearly all electrons are shock- accelerated to a power-law distribution of index p > 1, 
dn' e oc 7~ p d7 e , for j m < 7 e < jm, where n' e is the co-moving electron number density, j e is the electron 
random Lorentz factor, and we ignore a possible tail of thermal electrons with 7 e < j m . The maximum 
is determined by the synchrotron losses during the acceleration timescale, and, for most of the afterglow, 
is several orders of magnitude larger than j m . The minimum j m is set by parameterizing the total energy 
density stored in electrons after acceleration, as a fraction e e i of the internal energy density of the shocked 
fluid (which is given by the jump conditions at shock - Blandford & McKee 1976), and by the injection 
fraction £ of electrons that are accelerated at shock. The result is 



7n 



p - 2 1 - ( 7 M/7m) 1 " p 

P - 1 1 - (7M/7m) 2 " P 



1 



£ei ni p 
C m e 



(T-l) 



o>2 £el_ P - 2 m p 

C p - 1 m e 



(T-l) 



(10) 



m p and m e being the proton and electron masses, respectively. 

The comoving magnetic field B is assumed to be turbulent and is parameterized through the fraction 
£mag of the internal energy that is in the form of magnetic field energy, 

-l 1/2 



B = 



U 

'mSmag— — Bq 



jT+1 3 



-mao ^ 

7 



1 3a 3 -, 



(11) 
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where B = (87rr p d c 2 ) 1/2 = 1.94 ( i^o^o^) 1 ^ 2 G, no being the external medium particle density at the 
deceleration radius, in cm" - 3 , r , 2 = KT 2 r and the non-dimensional variables a — r/fd and u = U /M$c 
have been used. 

The radiative losses are given by a double integral over the remnant volume and the electron distribution, 

(dU) rad = - dV / dn' e { le )P' sy { le ) dt> , (12) 

J •'7m (f) 

which can be calculated for given B and dn' e ('y e ) at each point in the shocked structure. In equation (|l 
P'syile) — (4/3) aTh,c{B 2 /%tt)("{1 — 1) is the synchrotron power {axh being the cross-section for electron 
scattering) and dt' = {V 2 — l) -1 / 2 e _1 dr gives the comoving frame time. The electron distribution in each 
infinitesimal "sub-shell" within the volume of the shocked fluid is calculated by first initializing it at the 
time t' when the sub-shell is added to the shocked structure and then tracking the evolution of the electron 
Lorentz factor j e , subject to adiabatic and radiative losses: 

d 7e 1 dint" 1 a Th B 2 dt' 2 

T~ = ~Q~^ (7e-l)- 7 -Tile- 1 )- ( 13 ) 

da 3 da bir m e c da 



2.4. Differential Equations and Initial Conditions 

We can now add the contribution of the material injection at the reverse shock given by equation (g) 
to the evolution of T (eq. M ) , and obtain: 

dr _ iy[i - r 2 (i - p/3 f )]F e (r f ,T)u; M + (7 - i)«r M - 3(r 2 - ljiyV^ 

da a^ + ^r + ^^r + i)- 1 ^ ' 1 ' 

where fi = M/Mq is the dimensionless remnant mass and u>m = Mi n j/Mo is the dimensionless total injected 
mass. We can also obtain the differential equation for U by including the internal energy input (eq. ||) and 
the radiative losses (eq. Q): 

The first term in the numerator of equation (|l4|) and the first term in the right hand side of equation (|l5|) 
are switched off when Ty, calculated by integrating equation (Q), drops below r m . The differential equation 
for the mass of the remnant can be obtained with the aid of equations (^) and (||) : 

Solving the hydrodynamics of the remnant is therefore equivalent to integrating the set of coupled 
differential equations (Q), ([14]), (^5|) and (^6[). These equations are valid in any relativistic regime. If all the 
released ejecta has the same initial Lorentz factor, then the fireball is not entirely crossed by the the reverse 
shock if r < and the yet unshocked part of the ejecta and the shocked fluid move with different Lorentz 
factors. To avoid unnecessary complications, we simulate the dynamical evolution of the fireball starting 
from r = r<i and pass over the r < Td stage by making an appropriate choice of the initial conditions at 
r = Td- These initial conditions are determined by the definition of r^. The energy released before r^, can be 
safely neglected, so that, by equating the sum of the kinetic energy (1 + T^" 1 ) (Yd — IJMqC 2 and the lab frame 
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internal energy ~ r^r^ — 1)(Mq/Tq)c 2 at r d with the initial energy Eq = {Yq — 1)Mqc? , it is straightforward 
to show that T d = 0.62 T Q and U(r d ) = (T d - l)(M /r )c 2 ~ 0.62 M c 2 . Therefore the initial conditions are 

r(a = 1) =T f (a = 1) =0.62r , tt(a = 1) = 0.62, and /i(a = 1) = 1 + T 1 . (17) 



3. Analytical Asymptotic Light-Curves 

The temporal history of the afterglow flux received at Earth can be calculated analytically by assuming 
that the ejecta is either spherically symmetric or is a jet with axial symmetry, and that T is power-law in r. 
The last assumption is correct only over a certain range of times; a different treatment is needed when the 
remnant slows down to non-relativistic speeds. We consider here relativistically expanding remnants and, 
for simplicity, in this section we neglect energy injection and restrict our attention to the case when the 
remnant is adiabatic, as it is most likely that this stage lasts the longest (Waxman, Kulkarni & Frail 1998). 
Electrons can be either radiative or adiabatic. The former case is compatible with the assumption of an 
adiabatic remnant provided that electrons are not re-energized after shock acceleration or that e e z is small 
enough that most of the internal energy is stored in protons and magnetic fields and lost adiabatically. In 
what follows, we denote by 7 m the minimum Lorentz factor of the electrons (in the power-law distribution) 
that have just been accelerated, i.e. those electrons that are located very close the forward shock. 

For a relativistic and adiabatic remnant, to a good approximation, the Lorentz factor of the contact 
discontinuity evolves as (Meszaros et al. 1998) 

r ~ T d {r/r d y^l 2 , (18) 

where we consider only the a < 3 case. The definition of r d gives r d ex (Equ^T^ 2 ) 1 / 3, , where n d is the 
external medium particle density at r d , and where a multiplying factor that has a weak dependence on a 
has been ignored. For definitcness, we consider that the power-law behavior of the external medium density 
is manifested beyond a radius R d , up to which the external density is almost constant, with R d large enough 
to cover all the possible values of the deceleration radii encountered in fireballs with reasonable values of 
the parameters Tq and Eq. This approximation is not affecting the remnant evolution, as the afterglow 
radiation is emitted at radii much larger than R d . Using the relationship between the observer time T and 
the lab- frame time t, T oc t/T 2 , the T-dependence of the Lorentz factor is found to be 

T oc \(E /n d )^/ 2 r^] V(12 ~ 3Q) T -(3-a)/(8-aa) _ (19) 

Note that T is independent on Tq if the external medium is homogeneous. As seen by the observer, the 
(transverse) source size scales as TT if the ejecta is spherically symmetric. The received flux F v at the 
peak v p of the synchrotron spectrum is F Up oc (rT) 2 r 3 /^,, (Meszaros et al. 1998), where I' v , is the comoving 
synchrotron intensity at the comoving peak frequency v' p . 

Radiative Electrons. If electrons are radiative, then I' v , oc n' e (P' sy /v' p )t' sy oc n'^mB) -1 , where 
t' sy is the comoving synchrotron cooling timescale and where we substituted P' sv t' sv oc j m . The magnetic 
field can be calculated using equation (|ll|), where the comoving internal energy is determined by assuming 
that the lab frame internal energy of the adiabatic remnant is a (constant) fraction of the initial energy: 
Tu oc E a => B oc eUl g (n d - a T^ 2a E^ e r- a / 2 T. The same result can be obtained using the jump conditions 
at the forward shock: B 2 oc e mag (dU/dV') oc s mag p'T oc s mag p ex T 2 . By using equations jlC| ) (with p > 2), 
( |l9| ) and r oc T 2 T, one can calculate I' v , and the evolution of the observed peak frequency v v oc j^BT. If 
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observations are made at a frequency v < v p , the observer "sees" the low energy tail of the synchrotron 
spectrum, which has a slope of 1/3. Then F v — (y /v p ) 1 ^F Up , leading to 



F rad {v < v p ) oc e- m %X?t\^E^T . 
Above v p , the synchrotron spectrum has a slope —p/2, yielding: 

F rad {v > " P ) oce^-p^ e r 1 e- p Ei P+2 ^T-^-^ 



(20) 



(21) 



Note that F rad (v) is independent on the external medium parameters (a, rid) and on the fireball initial 
Lorentz factor and that it depends strongly (powers close to or above 1) on s e i and also on £ if v < v p and 
on Eq if v > v p . 

Adiabatic Electrons. If the electrons are adiabatic, then I' v , oc n' e (P' sy / v' p ) A' oc n' e BA', where A' is 
the comoving remnant thickness. The product n' e A' can be calculated using the fact that 4tt (n' e /()m p r 2 A' 
is the external medium mass swept up until radius r is reached. Below the spectral peak 



tadyv < Vp) oc e^ ag s el Q 



r -4a p 2(5-a) 2(3-a) 

1 o F Q n d 



1/(12— 3a) 



T 



(2-o)/(4-a) 



(22) 



For observations made above u p , the synchrotron spectrum has a slope — (p — l)/2, therefore 

1/(12-3q) 



-r-ia j 7 ,3(p+3)-a(3p+7)/4 2(3-a) 
1 h n d 



T -[3(p-l)/4]-[a/(8-2a)] _ ^3) 



Generally, the light-curve has a strong dependence on Eq and e e i, and also on £ if v < v p and on e mag 



if 



v > v p . Other dependences are weak to moderate. For a <; 2, the light-curve depends strongly on To and 
F a d{y < Vp) increases with T. The larger a, the faster F ad {v > v p ) decreases with T. 



Jets. If the ejecta is jet-like (Waxman et al. 1998), then equations ( |20D - (23) give the correct observed 
flux in the early afterglow, when the observer does not see the edge of the jet or the effect of the sideways 
escape of the ejecta (Rhoads 1998). For an observer located at an angle o b s relative to the jet axis and a jet 
of half-angular opening 9j et such that 6 b s -C 0j e t, the jet edge is seen after T drops below 9j et . In this case, 
the source size is oc r6j e t oc T 2 T9j e t- The light-curve of the afterglow from a jet-like remnant is given by: 



Tl a El- 2a r 



-(3- 



1/(12-3q) 



T l/(4-a) ^ 



r 2a p3(p+3)-a(3p+10)/4 -(3-a) 
L n d 



l/(12-3a) 



r -(3p/4)-[(2-a)/(8-2e0] 



F rad (u>v p )^s^-p/ i elr\ 2 - p 
if electrons are radiative and by 

TP d, „\r* cV3 c -2/3>5/3 f r -2a ^13-30^3-0] l/(12-3a) rp-l/U-a) 

i<a d \v < v p ) oc s^ ag e el C [1 n d J 

l/(12-3a) 



F ad (^>^)oc^+ 1 )/ 4 £ 7 1 C 2 ^ 



f -2q p 3(p+4)-«(3p+ll)/4 3-a 
1 h n d 



T -[3(p-l)/4]-[(6-a)/(8-2a)] 



(24) 

(25) 

(26) 
(27) 



if electrons are adiabatic. 



A comparison of equations (|2C|)-(|23|) and fl24|)-(|27D shows that the light-curve from beamed ejecta rises 
slower and decays faster than that from a spherical fireball. At the onset of the T < 6j^ t phase, the decay 
of the afterglow steepens by (3 — a)/ (4 — a), yielding a break in the light-curve. This phase lasts until the 
escape of the ejecta outside the cone in which it was initially released becomes important. Rhoads (1998) 
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has shown that in this case the remnant bulk Lorentz factor decreases exponentially with radius, and that 
the decay of the afterglow light-curve exhibits another break, but remains a power-law in the observer time. 
It can be shown that the time interval from the onset of the exponential phase and the beginning of the 
non-relativistic phase is 2.3 (9j et /10°) 2 times shorter than the duration of the T < 9~^ t phase and that unless 
Qjet 7° (which yields a very low probability of observing the afterglow) the sideways escape phase occurs 
after the remnant becomes non-relativistic. 

Mixed Electron Radiative Regimes. Equations (|2C|)-(|27|) were derived assuming that all the elec- 
trons are either radiative or adiabatic. The real situation is more complex, as the more energetic tail of the 
power-law distribution of electrons contains electrons that are radiative and contribute more to the received 
flux at some given frequency v 3> v p {"lm) than the less energetic 7 m -electrons, which become adiabatic early 
in the afterglow. In fact this is the case with most of the numerical X-ray and optical afterglows shown in 
the next section. If the 7,„-electrons are adiabatic, the flux at a frequency where the emission is dominated 

by more energetic and radiative electrons can be derived using the calculated for adiabatic electrons and 

p 

the fact that the spectrum has a slope — (p — l)/2 for frequencies above v v and below the peak frequency 
of the synchrotron emission from electrons that have a radiative timescale equal to the adiabatic one, and 
a slope — p/2 above this frequency. Interestingly, the result is the same as given by equations and ( p5| ) 
for radiative 7 m -electrons, i.e. only the constants of proportionality are altered. 

We should keep in mind that the above analytical derivations do not take into account the shape of the 
equal arrival time surface, i.e. the fact that photons that arrive simultaneously at detector were emitted at 
different lab-frame times. Moreover, we ignored the fact that there arc electrons with Lorentz factors below 
the 7,„ of the freshly accelerated electrons. For these reasons, the equations (^0])-(|27|) are of somewhat 
limited use and, for more accurate results, one must integrate numerically the afterglow emission. 



4. Numerical Results 

We have introduced so far the following model parameters: (1) dynamical parameters (Eg; rid, ct; Tq), (2) 
late energy injection parameters (Ei n j,T m , s) and (3) energy release parameters (£ mag ; SeliPi 0- To these 
one must add (9j e t,0 o bs) if the ejecta is jet-like. In this section we asses the effect of these parameters, 
and consider also the situation where Eq and Tq have an anisotropic distribution in the ejecta, which, in 
the simplest case, introduces one more parameter representing the angular scale of such anisotropy. We 
compare our numerical results to the observed X-ray (2-10 keV), optical (V magnitude) and radio (4.9 
GHz) afterglows. We will be looking in particular for the parameter values that yield X-ray and/or optical 
light-curves similar to GRB 970508, for which a fairly uniform time coverage is available. 



4.1. Spherically Symmetric Ejecta 

The simplest case is that of spherically symmetric ejecta with a single impulsive input of energy. Under 
the simplifying assumptions of a relativistic and adiabatic remnant, the equations ( |2C| ) ( p3| ) predict the 
asymptotic radio, optical and X-ray afterglow. For the range of times considered here, v p is below optical 
frequencies and only the radio emission shows a peak. This peak generally occurs before v v reaches few 
GHz, and it is due to the remnant's transition from the relativistic to the non-relativistic regime. For a 
homogeneous external medium (a = 0), radiative electrons and p = 2.5, the above-mentioned equations for 
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a relativistic remnant yield for v > v v (optical and X-ray fluxes) 

T? ^ =-1/8 = - 3 / 2 /— 1/2 7^9/8^-11/8 ( O o\ 

fo.x « s^ iag e e[ Q 1 h Q 1 1 , (28) 
while radio flux Fr is given by equation (p0|). If electrons are adiabatic, then 

Fo,x oc e^l/J^C^E^nTT^ , (29) 

Fn - elll^^ErnTT^ . (30) 

These analytical approximations are consistent, within their range of validity, with the numerical results 
shown in Figure 1. For all the the afterglows shown in Figure I, the 7 m -electrons become adiabatic for T 
between 0.01 and 10 days, while the remnant enters the non-relativistic phase at times between 10 and 300 
days, when a slow but steady steepening of the light-curves can be seen. Figure 1 also shows (with symbols) 
observational data taken from IAU Circulars, van Paradijs et al. (1997), Sahu et al. (1997), Frail (1997), 
Piro et al. (1998) or inferred from the data presented by Galama et al. (1997), Bartolini et al. (1998), and 
Sokolov et al. (1998). The numerical results are not meant to be fits to the observational data. 

If the optical and X-ray electrons are radiative, the afterglows arising from fireballs with larger initial 
energy or energy release parameters are brighter, as implied by equation (p8|). Fireballs with harder electron 
distributions lead to afterglows that have a shallower temporal decay (Figure 1 [a2] ) , as predicted by equations 
( pl| ) and (|23| ). If the electron injection fraction £ is sufficiently small, the radio afterglow can be undetectable 
(see eq. |3(|). The peak of the radio light-curve for the £ = 10~ 2 afterglow shown in Figure 1 is ~ 10 /iJy. For 
the same afterglow the synchrotron peak from 7 m -electrons remains above the optical range for several days, 
leading to an optical afterglow that is flat for the same duration (see Figure l[a2]). The non-detection of radio 
emission from a remnant that yields observable optical afterglows could also be due to an inhomogeneous 
external medium: the peak of the radio emission of the a = 2 (pre-ejected wind) case shown in Figure 1 is 
- 30 /zJy. 

The are some important differences between the light-curves arising from a fireball running into a 
homogeneous external medium and into a pre-ejected wind. First note that Figure l(bl) shows that when 
the electrons emitting at fixed frequency (here, in X-ray) are radiative, the afterglow is indeed independent 
of the external medium parameters n<j and a (if a < 1), as predicted by equation (pi]). The optical and 
the radio afterglows depend on a (this is also true for the X-ray light-curve if a > 1), indicating that in 
these cases the electrons that radiate most of the light in the corresponding energy bands are adiabatic (eqs. 
|P2| and [P3| ) . In a relativistic remnant, the lab-frame synchrotron cooling timescale t sy oc T/(^ e B 2 ) for 
electrons radiating at a peak frequency v p (^ e ) oc j 2 BT > v p (p/ m ) equal to a fixed observing frequency v is 
t sy oc y-^ 2 (T/B) 3 / 2 , leading to: 

t sy oc e^!XJ 2 ET'^- a)l \-^t^ . (31) 

which is constant in time for a homogeneous external medium, and increases as t z / 2 for a pre-ejected wind. 
The adiabatic cooling timescale increases as t, if the comoving density tracks the post-shock density. There- 
fore the electrons radiating above v v (i.e. in optical and X-ray) that are radiative, remain so during the entire 
afterglow if the external medium is homogeneous but eventually become adiabatic if the fireball interacts 
with a pre-ejected wind. The radiative regime of the electrons that emit at a given frequency changes with 
the index of the external medium, as implied by the increase of t sy with a (eq. [[n)), and as suggested by 
the light-curves shown in Figure l(bl) and l(b2): for a = the X-ray and optical electrons are radiative, 
for a — 1 only the electrons emitting in X-ray are radiative, while for a = 2 they are all adiabatic. Another 
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important difference between the homogeneous and pre-ejected external media models is manifested by the 
duration of the relativistic phase. From equation ( |l9| ) one can calculate the dependence on model parameters 
of the time T nr when the remnant becomes non-relativistic (T < 2): 

T nr oc (£ /^) 1/3 rif /(9 ' 3Q) . (32) 

Obviously, T nr is To-independent for a = 0, but it depends strongly on the fireball initial Lorentz factor in 

4/3 4/3 

the case of a pre-ejected wind: T nr oc T , implying that in this case the relativistic phase lasts <; 100 T 2 
times longer than in the homogeneous external medium case. The optical brightness of the a = 2 afterglow 
is correspondingly weaker, as shown in Figure l(b2). 

We have ignored the effects of low-frequency synchrotron self-absorption in the radio range, therefore 
the Figures l(a3) and l(b3) give essentially an upper limit to the optically thin radio flux expected in 
this case. A simple analytical derivation of the absorption frequency is straightforward (Meszaros & Rees 
1997), but it can easily lead to misleading results, since the fireball contains electrons with random Lorentz 
factors that span more than one order of magnitude, all emitting and absorbing the synchrotron radiation. 
Taking into account only the newly shocked electrons and ignoring a possible low-energy tail of the electron 
distribution below j m) it can be shown (Panaitescu & Meszaros 1998) that the self-absorption frequency is 
v a b ~ 6.4 (10 e mag) 1 / 5 (10 Eei)~ 1 nQ^ 5 Eq^ 2 T° GHz (at redshift 2 = 1) for a relativistic remnant and adiabatic 
electrons, where Jib, 52 = Eq/(10 52 ergs) and tiq — 71^/(1 cm -3 ). This result is valid until the remnant 
becomes non-relativistic or until the shocked material escape sideways, if the remnant is a jet. Therefore 
the optical thickness is r = 1.6 at 4.9 GHz for e e i = e mag — 0.1, uq = 1 and -E0.52 = 1, indicating that the 
radio fluxes shown in Figures l(a3) and l(b3) are overestimated by a factor of r(l — e~ T )~ 1 ~ 2. Post-shock 
mild re-acceleration of the cooling electrons or an electron (acceleration) injection fraction £ below unity can 
further decrease the radio flux by reducing the number of the low energy electrons in the remnant. 

In its simplest form considered in Figure 1, the fireball shock model obviously cannot explain departures 
from the power-law decay, such as observed in the optical afterglow of GRB 970508 near T ~ 2 days. A 
brightening of the afterglow may arise if there is a delayed energy input, as illustrated in Figure 2. The 
energy injection index s was set equal to a large value so that the input resembles a second relativistic shell 
that catches up with the initial fireball. For a delayed energy input Ei n j comparable to or larger than the 
energy of the remnant Eg, the light-curves exhibit a bump at the time of interaction between the two shells. 
The larger Ei n j is, the more prominent is the resulting bump. For lower T m , the collision takes place later, 
and this might explain a secondary departure from a power-law, apparent in the optical afterglow of GRB 
970508 at T Si 50 days, (the flattening of the light-curve could also be due to a constant contribution of the 
host galaxy - Pedersen et al. 1998). 

In Figure 2, the minimum Lorentz factor T m was chosen such that the numerical light-curve exhibits the 
brightening observed in the 970508 optical afterglow after T = 1 day. All light-curves shown in Figure 2 were 
calculated using the same fireball initial energy Eo = 6 x 10 51 ergs, delayed energy injection (from refreshed 
shocks) Ei n j = 3i?o (yielding a total energy Eo + Ej n j = 2.4 x 10 52 ergs), and T m = 11, and the same set 
of parameters (no, a; s mag ; £ e ;)- The model shown with dotted lines corresponds to constant parameters p 
and C, chosen such that the slope of the late optical power-law decay and the early time radio fluxes are 
close to the observed ones. The corresponding X-ray afterglow is too faint, while the early optical and late 
radio afterglows are too bright. Generally, such discrepancies cannot resolved by adjusting the dynamical 
parameters (Eo, Ei n j,T m ; no, a) or the energy release parameters (e mag ; £ e i), as changes in these parameters 
alter the multi-wavelength light-curves in a similar fashion. However, a physically plausible possibility is 
that changes occur in the parameters p and C which determine the shape of the synchrotron spectrum, and 
these can alter the light-curve in a given band without significant changes in other bands. 
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For times T <; 0.3 days in Figures 1 and 2, the synchrotron peak v p is below the optical band, so that 
the relative intensity of the optical and the X-ray fluxes is determined only by the slope of the spectrum 
above v p . This suggests that a brighter X-ray afterglow and a dimmer optical light-curve can be obtained 
by using a flatter electron index p, as illustrated by the early X-ray and optical fluxes shown with dashed 
lines in Figure 2. If p were held constant at 1.4 during the entire afterglow, the resulting optical light-curve 
would decay much slower than for p = 2.3 (see eq. plj), and thus would be clearly inconsistent with the 
observational data. A better simultaneous fit of the X-ray and optical afterglows can be obtained if one 
assumes that the electron index changes during the evolution of the remnant. In the model shown with 
dashed lines in Figure 2, we considered that the index p — 1.4 is constant until the second shell of ejecta 
catches up with the fireball (T ~ 2 days), and changes to p = 2.3 at the end of the collision between the two 
shells. The indices p before and after the delayed energy input were chosen so that the numerical result fits 
the early X-ray to optical emission ratio and the decay of the observed optical light-curve. The electrons that 
radiate most of the U-band light shown in Figure 2 are radiative, with some smaller contribution from the 
adiabatic 7 m -electrons, implying that the optical spectrum should have a slope close to — p/2. The change 
from p = 1.4 to p — 2.3 at T ~ 2 days is consistent with the optical spectral slopes reported by Djorgovski 
ct al. (1997) (-0.65 ± 0.30 at T ~ 1 day), Metzger et al. (1997) (-0.9 ± 0.3 at T ~ 2 days), and Sokolov et 
al. (1998) (-1.1 for T between 2 and 5 days). 

The radio afterglow at times shown in Figure 2 (T > 3 days) depends on the late value of index p. 
Unlike the emission at optical and X-ray energies, the emission at radio frequencies is due to all the electrons 
in the remnant, whether they are the first accelerated electrons (that have cooled and emit only in radio) 
or the more energetic, recently accelerated electrons (that radiate at higher frequencies but extend their 
emission down into radio through the low energy synchrotron tail of slope 1/3). The later electrons slightly 
dominate the radio emission after T ~ 10 days, and lead to the large fluxes shown with dotted line (£ = 0.2) 
in Figure 2. This contribution to the radio emission is diminished if the recently accelerated electrons have 
a higher post-shock acceleration Lorentz factor, which can be achieved if the electron injection fraction £ is 
decreased (see eq. |fL0| ). This is shown by the dot-dashed line in Figure 2, where it was assumed that the 
electron acceleration injection fraction drops from £ = 0.2 to £ — 0.05 when the remnant approaches the 
non-relativistic regime (r ~ 3). At the same time the optical afterglow exhibits a brightening due to the 
fact that for ( = 0.05 the synchrotron peak v v is closer to the optical range. 



4.2. Axially Symmetric Jets 

Jet-like outflows obviously reduce the energy requirements of fireballs, which, if extending over 47r sr, 
would require a total energy above 10 52 ergs to produce the optical fluxes observed in the afterglow of GRB 
970508. In Figure 3(a) we show light-curves arising from jet ejecta whose properties are isotropic within the 
opening angle 6j e t- From these numerical results, we can draw several conclusions: 

(1) As expected, the light-curve decay steepens when the observer sees the edge of the jet. This is shown 
by the departure of the dotted line (jet, observer located on the jet axis) from the thick solid line (isotropic 
fireball) around T = 6 days. The smaller 6j et , the earlier such a steepening occurs. 

(2) Jets seen at angle 9 b s < 9j e t do not exhibit the rise shown by jets with 9 b s > 9j et . 

(3) The larger o b s , the more delayed and dimmer the afterglow peak. For energies Eq 10 51 ergs, the 
optical emission from jets located at z = 1 that are seen at an angle 9 b s > 2 9j e t, is unlikely to be detected. 
The afterglow that fits best the observations is obtained when energy injection is included. The thin solid 
line in Figure 3(a) is for a total delayed energy input 4 times larger than the initial energy of the jet, leading 
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to a total available energy of 1.9 x 10 51 ergs. 

In a more realistic al scenario, the explosive event that generates the ejecta may lead to an angle- 
dependent energy distribution, as considered by Meszaros et al. (1998). Figure 3(b) (which is not meant as 
a fit to the afterglow of GRB 970508), shows the effect of such an anisotropic distribution for the particular 
choice where the energy per unit solid angle in the jet is an exponential in the polar angle 9: (dE n /dQ)(9) = 
(dEo/dQ) axis exp (—9/9e)- For 9e > the angular energy density decreases toward the jet edge while for 
9e < it increases. The same angular dependence (with the same angular scale 9e) was assumed for To. 
The initial Lorentz factor has no effect on the light-curve if the external medium is homogeneous, as shown 
in the previous section; the motivation for this choice was simply an isotropic mass distribution in the ejecta. 
To maximize the effect of the anisotropy in the ejecta, the observer was placed on the jet axis, and a large jet 
opening was chosen in order to separate this effect from the "edge effect" . In all cases, the energy density at 
9 = 0° was set to 10 52 / 71 " ergs/sr, which leads to the following total jet energies: E _ 52 = 1 for the isotropic 
distribution (9 E = oo), -E0.52 = 0.2 for 9 E = 9j et /3 and E ^ 2 = 8.6 for 9 E = —9j et /3. The light-curve decays 
agree qualitatively with the results of Meszaros et al. (1998): if dEo/dfl > 0, then more energy is emitted 
from fluid moving at larger angles relative to the LSC, arriving later at detector, and yielding shallower 
decays than in the isotropic case. Conversely, if dE /dQ < 0, then most energy is radiated away by the 
fluid moving close to the LSC; this radiation arrives earlier at the detector and leads to steeper light-curve 
decays. 

The case where the observer is located off the jet axis is considered in Figure 3(c). The parameters 
(dE n /d9) axis and 9e were chosen so that the total energy of the jet is the same in all cases. The conclusion 
that can be drawn from Figure 3(c) is that, for all other parameters fixed, the light-curve seen by an off-axis 
observer is determined mainly by the total energy of the jet and not by how this energy is distributed. 
The ironing out of the details of the angular energy distribution in an axially symmetric jet is due to the 
differential rclativistic beaming of the radiation emitted by fluid moving at angles between 9 b s — 9j et and 
9 bs + 9 jet relative to the LSC. 

Jets with the parameters given for Figure 3(a) and 3(c) can explain the rise and decay of the light-curve 
of GRB 970508 after T ~ 0.5 days. The emission detected in the early part (T 5s 0.5 day) of the optical 
afterglow may be due to some ejected material lying outside the main jet. In Figures 3(a) and 3(c) we 
show with dot-dashes lines the emission from such a large angle outflow, containing E n = 7.0 x 10 50 ergs, 
ejected isotropically outside of the central jet of opening angle 6j et — 10°, whose axis of symmetry is offset 
by 9 bs = 14° relative to the LSC. The sum of the light-curves from such a two-component ejecta (central 
jet and large angle outflow) matches well the features observed in the afterglow of GRB 970508. The X-ray 
afterglow can be fitted as before together with the optical, by making an appropriate choice of the electron 
index p in the jet and in the large angle outflow. 

5. Conclusions 

Previous models of GRB afterglow light-curves from cosmological fireball shocks (e.g. Meszaros et 
al. 1998; Rees & Meszaros 1998; Sari, Piran & Narayan 1998) have used analytical descriptions based 
on scaling laws valid in the asymptotic limits. These require simplifying assumptions and involve various 
undetermined parameters. The most important analytical results on the afterglow light-curve are given in 
section §3. They should be used with care when making comparisons with observed power-law decays, as 
electrons with different random Lorentz factors can be in different radiating regimes. Generally, those elec- 
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trons radiating in optical and X-ray are radiative, while those radiating at radio frequencies are adiabatic, 
at least as long as the remnant is relativistic. Moreover, the analytical light-curves do not take into account 
the shape of the equal-arrival time surface, and assume that there is a one-to-one correspondence between 
the lab-frame time of emission and the detector time. 

Numerical calculations provide the environment where the assumptions made in analytical derivations 
can be tested and relaxed, and where results are expected to be more accurate. In some cases, like that of a 
fireball in a mildly relativistic regime, or like that of jet ejecta seen at an angle o b s ^ 0°, it is cumbersome 
to obtain analytical results. At the level of numerical calculations, effects arising from the viewing geometry 
(the equal arrival time surface is not the same as the equal lab-frame time surface) or from details of the 
energy release (e.g. an accurate tracking of the evolution of the electron random Lorentz factor j e ) can be 
properly accounted for. 

We have solved the differential equations for the afterglow evolution and integrated the remnant emission 
to calculate light-curves with different model parameters. Energy injection (refreshed shocks), angular 
anisotropy and jet-like structure of the ejecta allow for a variety of possible behaviors of the numerical 
light-curves, even under the assumption of axial symmetry in the remnant. More than one scenario could 
explain a fairly large fraction of the optical data of the GRB 970508 afterglow. A spherically symmetric 
ejecta with energy injection up to a total energy of 2.4 x 10 52 ergs, or a jet of opening 10° seen at an angle 
of 14°, in which energy is injected up to a total of ~ 2 x 10 51 ergs, both located at redshift z = 1, seem to fit 
most the mentioned afterglow. Such energies arc quite conservative in a cosmological scenario, and clearly 
do not require any drastic departures from the simple fireball/firejet scenario. Using a variable index of the 
electron power-law distribution, we obtained a simultaneous good fit of the X-ray and optical afterglow of 
GRB 970508. Synchrotron self-absorption, post-shock re-energization of the electrons or a decrease in the 
electron acceleration fraction reduce the radio fluxes obtained numerically and yield a better fit of the radio 
data, but this comparison is much more uncertain and model-dependent. 

The first generation of fireball shock models of afterglows were characterized by great simplicity and 
have predicted power-law decaying light-curves. As one would expect, relaxing some of the assumptions that 
are usually made in the simplest versions of these models, such as isotropy of the ejecta or constancy of the 
parameters that quantify the energy release, leads to an improved agreement between numerical results and 
observations. All of the models presented here still contain simplifying assumptions (e.g. axial symmetry, 
power-law delayed energy input), which were taken as a starting point in investigating the features of the 
numerical light-curves. While the present data do not require it, relaxing these assumptions could lead to 
even more diverse afterglow light-curves. The variety of behavior exemplified by the models we have discussed 
highlights the potential importance of afterglow data as diagnostics for the dynamics and anisotropy of the 
ejecta, and emphasizes how much more can be learned when the sample has grown larger. 

We are grateful to NASA NAG5-2857, NAG 3801 and the Royal Society for support, and to R.A.M.J. 
Wijers and D. Reichart for useful comments. 
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Fig. 1. — The effect of the energy release (left graphs) and dynamical (right graphs) parameters on the light- 
curve from a simple spherically symmetric fireball, without refreshed shocks. Legends give the parameters 
for each case. For each curve there is another one that differs in only one parameter, allowing assessment of 
its effect. Other parameters are: i?o,52 = 1, no = 1, a = for graphs (al)-(a3), and e ma g = 1CP 2 , e e i = 0.1, 
p = 2.5, ( = 1 for graphs (bl)-(b3). Observational data: open symbols are for GRB 970228, filled symbols 
for GRB 970508. Graphs (a2) and (b2): V magnitudes inferred from Rc magnitudes are shown as squares. 
Error bars are given only for magnitude errors larger than 0.5. Graphs (a3) and (b3): triangles indicate 
upper limits. The radio light-curves for the £ = 10~ 2 and a — 2 afterglows have peak fluxes of 10 /xJy and 30 
yuJy, respectively, and do not appear in graphs (a3) and (b3). A redshift z = 1 in a H = 75 kms _1 Mpc _1 , 
= 1 Universe is assumed. The radio fluxes plotted are the optically thin upper limits; the inclusion of 
synchrotron self-absorption and/or electron rc-encrgization would lead to lower radio fluxes (§4). 
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Fig. 2. — Effect of refreshed shocks in an isotropic fireball, caused by a late energy input which is a power 
law in the Lorentz factor Tf of the ejecta that catches up with the fireball (§2.2). All models have the same 
initial and injected energies £'0,52 = 0.6, Ei n j = 3Eq, as well as the same minimum Lorentz factor of the 
delayed energy input T m = 11. The injection index s has a large value, leading to an impulsive energy 
input at r m and to a distinctive step- like brightening of the afterglow. Other parameters are: e mag = 0.1, 
e e i = 0.1, no = 1, a = 0, z = 1. An absorption of Ay = 0.25 mag (Reichart 1998) at the source redshift 
was assumed. The electron index p and acceleration fraction ( for each model are given in the legend of the 
optical light-curves. They are constant for the model shown with dotted lines, p changes at the end of the 
delayed energy input for the dashed and dot-dashed lines models, while £ decreases when the remnant ends 
the relativistic expansion only for the model shown with a dot-dashed line. Symbols represent the data for 
the GRB 970508 afterglow. 
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Fig. 3. — Optical light-curves from jet-like ejecta, compared to data points for GRB 970508. (a) An 
outflow which is isotropic within a jet of opening half-angle 9j et — 10°, seen at the different angles 9 a b Sl 
for E = 3.8 x 10 50 ergs, n = 1, a = 0, e mag = 0.1, e e i = 0.1, p = 2.5, ( = 1, z = 1. For comparison, 
the afterglow from a spherically symmetric remnant with the same parameters, except Eq = 5 x 10 52 ergs 
(yielding the same energy density per solid angle), is also shown (solid thick line). A numerical light-curve 
matching the observational data (solid thin line) corresponds to o b. s = 14° and energy injection characterized 
by Ei n j — 1.5 x 10 51 ergs, T m — 2 and s = 1.5. (b) Effect of an anisotropic angular distribution of energy 
inside a jet with 0j et = 60°, 9 i, s = 0°, (dE n /dQ) axis = l0 52 /ir ergs/sr. Other parameters (n 07 a;e mag ,e e i 7 p) 
are the same as for graph (a). The legend gives the angular scale 9e (see text), (c) The same jet as in (a) 
seen at 9 b s = 14°, but with different energy per solid angle distributions. All jets have the same energy 
Eq = 1.5 x 10 51 ergs, isotropically distributed (solid line), exponentially decreasing toward the jet edge 
(dotted line) or exponentially increasing toward the edge (dashed line). Also shown in graphs (a) and (c) 
with dot-dashed lines is the contribution from an ejecta which is isotropic everywhere outside of the jet with 
opening angle 8j et = 10° and orientation 9 b s = 14°. The isotropic component has an energy 7 x 10 50 ergs 
(other parameters are as for [a]) and can account for the early (T ^ 1 day) afterglow emission. 



